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We consider the effect of tlie violation of tfie equivalence principle (VEP) by the massive neutrino 
component on the Cosmic Microwave Background angular power specrum. We show that in the 
presence of adiabatic and isocurvature primordial density perturbations the Planck surveyor can 
place limits on the maximal VEP by the massive neutrino component at the level of few x 10^^, 
valid in the general relativity, for the case in which the gravity is the single source of VEP. This 
work has been performed within the framework of the Planck/LFI activities. 
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The atmospheric neutrino experiments [Q, provide 
strong evidence for neutrino oscillations implying a non- 
zero neutrino mass with a lower limit ~ 0.04 — 0.08 eV. 
The latest measurements of the contribution of double 
beta decay to the neutrino mass matrix |^ place an up- 
per limit on the neutrino mass of < 0.26eV. The 
direct implication of these results is the non-negligible 
contribution of the hot dark matter (HDM) to the to- 
tal mass density of the universe and the existence of 
three massive neutrino flavors (i.e. a density parame- 
ter n^h'^ « ^^^^ m,/93eV, where h = Hq/WO Km s^^ 
Mpc~^ is the dimensionless Hubble constant). The in- 
troduction of a HDM component in the form of mas- 
sive neutrinos with the total mass in eV range is also 
required for consistency between the Cosmic Microwave 
Background (CMB) anisotropy at small scales with the 
Large Scale Structure (LLS) of the universe derived from 
galaxy surveys [|, |, |, §. 

The most popular mechanisms examined as possible 
solutions to the neutrino oscillation problem are the 
mass-mixing vacuum and the matter induced (MSW) os- 
cillations |8|, having as basic underlying hypothesis the 
existence of non-degenerated neutrino masses and the 
non equivalence of neutrino weak and mass eigenstates. 

An alternative mechanism to these approaches is the 
neutrino oscillations due to the Violation of the Equiva- 
lence Principle (VEP) that can occur if the week equiv- 
alence principle is not satisfied, i.e. the coupling of neu- 
trinos with the gravitational field is non-universal and 
the flavor eigenstates are not identical to the states that 
couple to gravity ^ Such mechanism has 

been investigated by using the experimental results from 
Super-Kamiokande ( SK) [ p^ Q and Sudbury Neutrino 
Observatory (SNO) [|5) collaborations. Recent works 
have considered the possibility to test VEP by using as- 
trophysical neutrino sources as well as the effect of 
neutrino oscillations induced by VEP as an explanation 
of the pulsar kick velocity [jl^ . 

A deviation from the equivalence principle can be 
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parametrized by assuming that the parametrized-post 
Newtonian (PPN) parameters for a given metric have 
particle-dependent values In the conformal New- 

tonian gauge the perturbations in the Robertson- Walker 
spacetime are characterized by two scalar potentials ^' 
and $ which appear in the line element as p9[ : 

ds^ ^ a^(T)[-(l + 2a^')dT2 + (1 - 2-/<^)dx'dx^] , (1) 

where a(T) is the cosmic scale factor (ao=l today) and 
T is the conformal time, a and 7 are theory dependent 
parameters: in the general relativity a = 7 = 1 while 
in alternative theories of relativity satisfying the equiv- 
alence principle a and 7 need not to be unity but their 
values are universal (the same for all kinds of particles). 
In eq. (1) ^' plays the role of the gravitational poten- 
tial in the Newtonian limit Limits on the VEP for 
a = 1 have been obtained from the supernova SN1987A 
by comparing the arrival time for photons and neutrinos 
(17^—7^1 < few xlO^^) and for neutrinos and 

antineutrinos ([7^^ — 7pJ < 10 ) p^, |25[|. 

The CMB anisotropy holds the key for understanding 
the seeds of the cosmological structures and allows the 
measurement of the most important cosmological param- 
eters. The new generation of CMB experiments such as 
MAP Q and Planck Q will achieve enough precision 
to reveal the structure formation process up to arcminute 
angular scales. The aim of this paper is to study the 
possibility to test with Planck the VEP by the massive 
neutrino component. 

The evolution of the density fluctuations from the 
early universe involves the integration of coupled and lin- 
earized Boltzmann, Einstein and fluid equations plf for 
all the relevant species (e.g., photons, baryons, cold dark 
matter, massless and massive neutrinos). The differences 
introduced in the gravitational potential due to the VEP 
by the massive neutrino component can generate met- 
ric perturbations that affect the evolution of the density 
fluctuations of all the components, leaving imprints on 
the CMB angular power spectrum. 

A recent work analysed the possibility to detect 
neutrino oscillations with accurate CMB experiments. In 
the presence of a large neutrino asymmetry the mass- 
mixing vacuum neutrino oscillations left detectable im- 
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prints on the CMB anisotropy and polarization power 
spectra. For a neutrino flavor Vi with energy Ei,^ and 
mass m^. < 0.26 eV, the survival probability after the 
propagation through the distance L in the expanding 
universe, V^^^^^-iE^., L) = 1 - sin^ 26*0 sin^(7ri/A), 6*0 be- 
ing the vacuum mixing angle, does not exhibit a depen- 
dence on the oscillation length A ^ E^.^, as L >> A and 
< sin2(L/A) >^ 1/2. 

For the VEP induced oscillation, the oscillation prob- 
ability is given by the same equation, replacing 9q with 
the gravitational mixing angle, 9g, and with the oscilla- 
tion length given by A ~ {'^Aa^-Ei,.)~^ , where 5' is the 
gravitational potential and Aa^- — I ~ ai parametrizes 
the VEP by the neutrino flavor [|l2|, Also in the 
VEP case, L >> A for neutrinos with masses < 0.26 
eV and the oscillation probability is not sensitive to the 
VEP parameter Aa^. It follows that the mass- mixing 
vacuum oscillations can not be distinguished from VEP 
induced oscillations by using the CMB anisotropy. 

On the other hand, the imprint of the clusterization 
process on the CMB anisotropy depends on the dynam- 
ics of the involved particles ||2^ and can be then sen- 
sitive to the VEP. We evaluate the effect of the VEP 
by the massive neutrinos computing the CMB temper- 
ature fluctuations generated in the non-linear stages of 
the universe evolution when clusters and superclusters of 
galaxies start to form. We place limits on the maximal 
VEP by the massive neutrinos in the general relativity, 
valid for the case in which gravity is the only source of 
VEP. 

In we studied the CMB anisotropy induced by 
the non-linear perturbations in the massive neutrino den- 
sity associated with the non-linear gravitational cluster- 
ing process. This process generates metric perturba- 
tions leading to a decrease in the CMB anisotropy power 
spectrum of amplitude AT/T « 10^^ for angular reso- 
lutions between ~ 4 and 20 arcminutes, depending on 
the cluster mass scale Mdust and the neutrino fraction 
fy = f7,y/(rib + 51crfm), where D.^-, and 0.cdm are the en- 
ergy density parameters for neutrinos, baryons and cold 
dark matter particles, respectively. In the Newtonian 
limit, the neutrino gravitational clustering can be de- 
scribed as a deviation from the background by a potential 
^ given by the Poisson equation: 

V2*(f,a) = 47rG'aVm(a)^m(r,a) , (2) 

where f is the position 3- vector, is the matter den- 

sity and 6m(f,a) is the matter density fluctuation. The 
equations governing the motion of each particle species 
i (i = cold dark matter, baryons, neutrinos) in the ex- 
panding universe are given by |2^, |2^ : 

^ = -ai7(a)a,W, ^ ^ q,^[a^ H{a)]-\ (3) 
da da 

where qi is the comoving momentum and H[a) is the 
Hubble expansion rate: 

H^{a) = ^^{flm/a^ + Qr/a'^ + J^a + ^k/a^) ■ 



Here G is the gravitational constant, flm = + ^cdm + 
fli, is the matter energy density parameter, fir, and 
are the density parameters for radiation 
(including the contribution from photons and relativistic 
neutrinos), vacuum energy and the curvature of the uni- 
verse. The parameter ai from the eq. (3) accounts for 
the coupling of each z*'' particle species with the gravi- 
tational potential. The Newtonian description given by 
eqs. (2) and (3) applies in the limit of the weak gravita- 
tional field if, at each time step, the size of the non-linear 
structures is much smaller than the causal horizon size 
(i.e., the background curvature is negligible). 

By using numerical simulations [20| , we evaluate 
the effect of the gravitational clustering on the CMB 
anisotropy power spectrum first assuming ai = I for all 
the particles and then considering that the equivalence 
principle is violated only by massive neutrinos (a^ 7^ !)■ 

We consider a flat ACHDM cosmological model with 
f]„=0.38, f}A=0.62, Ho=62 Km s"^ Mpc"\ three mas- 
sive neutrino flavors giving a neutrino fraction /y=0.06 
(TOi/=0.78eV, ri^=0.022), and adopt a cluster mass 
Mciust =1.5xlO^^M0. We assume a primordial power 
spectrum with adiabatic perturbations and a scalar spec- 
tral index n'"'' = 0.98 and neglect, for simplicity, the con- 
tributions from the tensorial modes (gravitational waves) 
and the reionization effects. This model is consistent with 
the LSS data and the latest CMB anisotropy measure- 
ments, allowing at the same time a pattern of neutrino 
masses consistent with the results from neutrino oscilla- 
tion and double beta decay experiments. Panel a) from 
Fig. 1 presents the CMB anisotropy power spectra ob- 
tained from numerical simulations for different values of 
Aui, = \ — a^. The VEP parameter Aui, ^ induced 
metric perturbations, leading to a decrease of the CMB 
temperature anisotropy at small scales as Aoii, increases. 

In order to evaluate to what extent the Planck sur- 
veyor can test the VEP by the massive neutrino compo- 
nent, we compute the 1 — a errors on the estimates of 
the cosmological parameters and on the parameter Aun 
by using the Fisher information matrix approach. We 
consider for this computation only the Planck "cosmo- 
logical" channel between 70 and 217 GHz, a sky coverage 
fsky = 0.8 and neglect for simplicity the foreground con- 
tamination [ p9[ . 

The CMB anisotropics are sensitive to the cosmologi- 
cal parameters through the evolution of the cosmological 
perturbations since the end of the inflation. The sim- 
plest one-field inflation model predicts primordial adia- 
batic density fluctuations. Isocurvature density pertur- 
bations can also arise in multiple field inflation mod- 
els |3^, ^ IUJ. In order to test the sensitivity of 
the CMB anisotropy to the VEP parameter Aa^, for 
conditions more general than the pure adiabatic case, 
we include in our Fisher analysis the parameters for 
the amplitude of the baryon isocurvature density mode 
< lb h >, neutrino isocurvature density mode < 1^ Iv >, 
their cross-correlation < lb > as well as their cross- 
correlations with the adibatic mode < Alb > and < 
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FIG. 1: Panel a): the CMB anisotropy power spectra (only adiabatic density mode) including the non-linear gravitational 
clustering. Different values of Aqi, = 1 — are considered: Aa^ — (solid line), Aa^/ = 2.5 x 10~^ (dots), 2.5 x 10~^ (dash- 
dots), 1.8x 10~^ (dash-three dots), 1 x 10~^ (long dashes). We report also the case in which the non-linear gravitational clustering 
and the VEP are not included (dashed line). Panel b): relative differences, in terms of [1 — C,^''^(Aa^ > 0)]/C,^''^(AQf^ = 0), 
between the case Aq^ = and the cases with Aa^/ > identified by the same lines as in panel a). Panel c): the CMB 
anisotropy power spectra of the four considered fiducial models compared to recent CMB anisotropy data: FMi (solid line), 
FM2 (dashes), FM3 (dash-dots), and FM4 (dash-three dots). Power spectra normalized to COHE-DMR 4- year data j30| (see 
also the text). 



TABLE I: 1 — (T errors on the cosmological parameters, ampli- 
tudes of different modes and VEP parameter Aq^ potentially 
measurable with Planck (see also the text). 
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AIu defining the following parametric fiducial power 
spectrum | |35| ]: 

CM = Cf{p) (4) 
+ C/" (p)+ <hh> Cl' (p) 

+ <AI,> Cf^'"{p)+ <Ah> Cf'ip) 
+ </,/b>C/'-'^(p), 

where p represents the cosmological parameters of our 
ACHDM cosmological model, Cf{p}, Cf-ip), Cl''{p) are 
the power spectra for the adiabatic, neutrino isocurvature 



and baryon isocurvature density mode, C/*'^" and C/*'^'' 
are the power spectra of the cross-correlation between 
the adiabatic and neutrino isocurvature density mode 
and adiabatic and baryon isocurvature density mode and 
Cl"'^''{p) corresponds to the cross-correlation between 
neutrino isocurvature and baryon isocurvature density 
mode. The power spectra Cf^, Cj"'^" and C/'"^'' are 
computed by using the CMBFAST code with pri- 
mordial conditions including also the neutrino isocurva- 
ture density mode given by p^ . The cross-correlation 
power spectra are computed by running the CMBFAST 
code with two modes excited and then substracting the 
power spectra computed when each mode was individu- 
ally excited, following the prescriptions of js^. For each 
mode we assume a primordial power spectrum with a 
scalar spectral index n""** = n^^° — n]^° — 0.98. For 
this computation we do not consider the cold dark mat- 
ter isocurvature density mode, as its power spectrum is 
very similar to that obtained for the baryon isocurvature 
density mode. We consider four fiducial models: only the 
adiabatic mode (FMi), the adiabatic mode, the neutrino 
isocurvature density mode and their cross-correlation 
(FM2), the adiabatic mode, the baryon isocurvature den- 
sity mode and their cross-correlation (FM3), and the 
adiabatic and all isocurvature modes as well as their 
cross-correlations (FM4). Panel c) of Fig. 1 presents the 
four fiducial power spectra obtained for < /;,/;, >= 0.1, 

< >= 0.1, < Ah >= 0.07, < Ah >= 0.07, 

< /i/ /(, >= 0.05 and assuming Aa^=0: they turn to be 
in a rough agreement with the current CMB anisotropy 
data [37-40]. 

Table 1 presents the 1 — ct errors on the cosmo- 
logical parameters, amplitudes of different modes and 
the VEP parameter Aa^ potentially measurable by us- 
ing the Planck observations of the CMB temperature 
anisotropy. Although the Fisher matrix approach is sen- 
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sitive to the considered fiducial model and to the number Romania Bilateral Agreement (Project n. 37) financial 

of parameters involved in the computation, this results support. 

shows that the Planck surveyor will be able to place an 

upper limit on the VEP parameter Aajy of few x 10^*^, 

valid in the general relativity, for the case in which the 

gravity is the only source of VEP. 
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